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TECHNICAL MEMORANDUM X-53953 


A SIMPLIFIED DEAD RECKONING NAVIGATION SYSTEM 
FOR THE MANNED LUNAR ROVING VEHICLE 

INTRODUCTION 


Problems related to lunar surface exploration have been studied by 
MSFC for the past several years. Early in these studies, a vehicle capable 
of navigating over the lunar surface was required. Subsequently, several 
industrial firms proposed and fabricated candidate lunar surface vehicles. 

In addition, analytical studies concerning lunar navigation and exploration have 
been carried out [1-6], 

More recently, plans are being developed for exploration in the 1972- 
1980 time-frame. Current work is being directed to the development of two 
vehicles. The first one, used in early exploration, is the manned lunar roving 
vehicle (MLRV) . This system will be driven by an astronaut onboard and will 
require dead reckoning navigation. The second system, used in advanced 
exploration, is the dual mode lunar roving vehicle (DMLRV) . This system 
will be driven by either an astronaut onboard or by earth-controlled remote 
methods. In addition to dead reckoning navigation, a position-fix system will 
also be required on the DMLRV, 

Many concepts for dead reckoning navigation have been proposed [5] . 
The navigation hardware for implementing these concepts includes various 
combinations of vertical gyros, directional gyros, pendulous inclinometers, 
accelerometers, sun sensors, star trackers, doppler radar, IR and rf earth 
trackers, and odometers. 

This paper gives a technique for dead reckoning navigation of an MLRV. 
The hardware involved consists of two free gyros hard- mounted to the frame 
of the vehicle, odometers, and tachometers. For the system to be self- 
contained, some form of signal processing (small computer) will also be 
needed. 



The salient requirements of the MLRV navigation system are: 


1. Navigation within a 5-km radius of the lunar excursion module (LEM) . 
Within the 5-km radius the lunar roving vehicle (LRV) will travel up to 30 km. 

2. An indication of the total distance traveled, with an accuracy of 
±2 percent, must be available. 

3. At a 5-km radial distance, the heading to the LEM must be known 
within ±3 degrees. 

4. The navigation system must be capable of returning the LRV within 
500 meters of the LEM. 

5. The system must have the capability of revisiting an established 
site. (Without visual cues this may be difficult to accomplish since an 
accuracy of better than 500 meters would be required. ) 

6. The LRV must be capable of sorties of a maximum 3-hour duration. 

7. The LRV must be capable of a 3-day period of lunar exploration. 


CHARACTERISTICS OF THE PROPOSED NAVIGATION SYSTEM 


First it should be stated that there are many navigation schemes, with 
better accuracy than the one proposed, which will satisfy the above require- 
ments. However, high resolution navigation systems carry with them the 
penalty of an increase in weight, volume, and complexity. The proposed two- 
free-gyro system is easy to implement and relatively light in weight. 

The system is defined by the following characteristics. At time = 0 
(the start of any given exploration trip) , vehicle body axes and moon-fixed 
inertial axes are aligned with one another. This can be accomplished by caging 
the gyros and bringing the gyro rotors up to full speed while in the caged 
condition. The orientation of the coordinates is arbitrarily defined by the 
alignment of the vehicle with any outstanding landmark. It is desirable that 
this landmark be recognizable from local terrain photographs. The center of 
the coordinate systems is on the surface of the moon and in the general vicinity 
of the shelter (LEM) site. As the moon rotates around its polar axis, the 
origin of the moon-fixed inertial coordinate system also rotates with the surface 
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of the moon. The total velocity of the vehicle, relative to the surface of the 
moon, is resolved into components along the X, Y, and Z axes of the moon- 
fixed inertial system. These components are integrated to give the distance 
traveled relative to the moon coordinate system. For the system to be com- 
pletely accurate the moon-fixed reference system should be Schuler-tuned to 
remain fixed relative to the moon's surface. However, during any 3-hour 
sortie the moon will rotate a maximum of approximately 1. 5 degrees. 
Neglecting this rotation will introduce small errors in the reference bearing 
of the vehicle, the. heading of the vehicle with respect to the shelter area, 
and the distance traveled relative to the moon reference system. As stated 
previously, Schuler tuning can be employed to eliminate these small errors, 
but the added complexity is not justified. The MLRV navigation requirements 
can be met using state-of-the-art gyros (with their associated drift rates) 
without Schuler tuning. 


COORDINATE SYSTEM DEFINITION 

Assume that the vehicle body coordinate system is defined so that the 

X axis is along the aft-forward sense with X positive in the forward direction. 

The Z axis is perpendicular to the X axis and positive upward. The Y axis 

completes a right-hand orthogonal coordinate system. This set of axes will 

be denoted by the subscript "B" to indicate body axes (X , Y , Z ). A 

B B B 

sketch of the coordinate system located in the vehicle frame is shown in 
Figure i. 



T B 


Figure 1. Vehicle body coordinate system. 
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An inertial coordinate system with its center at the LEM site and 
whose origin moves with the moon rotation is defined so that it is initially 
aligned to the body coordinates at t = 0. The coordinates of this system are 
given by Xj, and Z . 

At any time, the body coordinate system can be located with respect to 
the Xj, Yj, Zj. system by a set of sequential rotations. Many such sets of 

rotations exist that can serve this purpose. However, by selecting the 
sequential rotations in a certain way, information can be obtained which will 
allow one to determine the body axes location (relative to the X^, Y^, Z 

frame) by direct readings of gimbal pickoffs from the gyros. A method for 
accomplishing this is given later. 

In the work that follows, the mathematical equations for the navigation 
system are derived. The variables for these equations are related directly to 
the gimbal angles of the gyros. Thus knowledge of the vehicle velocity and the 
gimbal angles are sufficient for determining the vehicle position, bearing, and 
heading with respect to the shelter site. 


DEVELOPMENT OF THE NAVIGATION SYSTEM EQUATIONS 

Consider the schematic of two free gyros mounted to the vehicle frame 

as shown in Figure 2. Initially the X , Y , Z frame is aligned to the 

B B B 

X , Yj, Zj. frame. Four gimbal angles are defined on the diagram. Two 

transformations relating the X , Y , Z axes to the X Y Z axes are 
« j* i *i 13 13 111 

given as follows. 


X 


X 

Y 

= [A2(Z) ] [Ai (Y) ] [U (X) ] 

Y 

Z 

B 

Z 


4 
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and 



[A2 (Z) ] [A1 (Y) ] [U (X) ] = [A4 (X) ] [A3(Y)] [L(Z)] (4) 


The transformation components 

are given 

by 





cos(Al) 

0 

-sin(Al) 


cos(A2) 

sin(A2) 

0 

A1(Y) = 

0 

i 

0 

;A2(Z) = 

-sin(A2) 

cos (A2) 

0 


sin(Ai) 

0 

cos(Al) 


0 

0 

1 



cos (A3) 

0 

-sin (A3) 

i 

0 

0 

A3 (Y) = 

0 

1 

0 

; A4(X) = 0 

cos (A4) 

sin(A4) 


sin (A3) 

0 

cos (A3) 

0 

-sin(A4) 

cos (A4) 


6 




1 

0 

0 

cos(L) 

sin(L) 

0 

U(X) = 

0 

cos(U) 

sin(U) ; L(Z) = 

-sin(L) 

cos(L) 

0 


0 

-sin(U) 

cos (U) 

0 

0 

1 


For simplicity, the notation sin(Ai) = Si, cos(Ai) = Cl, sin(U) = SU, etc. 
is adopted. 

Note that the gimbal angle readings Al, A2, A3, and A4 are actual 
sequential angles. Angles U and L cannot be read directly from the gimbal 
angles but can be calculated (as shown later) in terms of Al, A2, A3, and 
A4. 

From equation (4) the following expression may be written. 

[U (X) ] [L(Z)]" 1 = [Al(Y)]- 1 [A2(Z)]-* [A4(X)] [A3(Y)j (5) 



Cl 

0 

SI 

C2 

-S2 

0 

1 

0 

0 


C3 

0 

-S3 

[U(X)][L(Z)]“ 1 = 

0 

1 

0 

S2 

C2 

0 

0 

C4 

S4 


0 

1 

0 


-SI 

0 

Cl 

0 

0 

1 

-0 

-S4 

C4 


S3 

0 

C3 


C1C2 

-C1S2 

SI 


C3 

0 

-S3 

S2 

C2 

0 


S4S3 

C4 

S4C3 

-S1C2 

S1S2 

Cl 


C4S3 

-S4 

C4C3 
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giving 


iU(X)][l<(Z)] 


Also, 


-l 


C 1C2C3 - C 1S2S4S3 + S1C4S3 ~C 1S2C4 - S1S4, 

S2C3 + C2S4S3 C2C4 


(-S1C2C3 + S1S2S3S4 + C1C4S3 

S1S2C4 

- C1S4 


CL 

-SL 

0 

[U(X)][L(Z)]~ 1 = 

CUSL 

CUCL 

SU 


-SUSL 

-SUCL 

CU 


-C1C2S3 - C1S2S4C3 + S1C4C3 
-S2S3 + C2S4C3 
S1C2S3 + S1S2S4C3 + C1C4C3 


( 6 ) 


(7) 


Therefore, equating (6) and (7) yields 

CL = C1C2C3 - C1S2S4S3 + S1C4S3 
SL = C1S2C4 + S1S4 
SU = -S2S3 + C2S4C3 
CU - S1C2S3 + S1S2S4C3 + C1C4C3 


( 8 ) 

(9) 

( 10 ) 

(ID 


Let the sequence of rotation be given by equation (2) . Then 


X 


C3C1 

C3SL 

-S3 

Y 

= 

-C4SL+ S4S3CL 

C4CL + S4S3SL 

S4C3 

Z 

B 

S4SL + C4S3CL 

-S4CL + C4S3SL 

C4C3 


X 

Y 
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This sequence of rotations is shown in Figure 3. 



Figure 3. Sequential angle rotation defining the vehicle body axes 
with respect to the inertial axes. 

The Xj, Yj, Zj components are related to the body components by 

X C3CL -C4SL + S4S3CL S4SL + C4S3CL X 

Y = C3SL C4CL + S4S3SL -S4CL + C4S3SL Y 

Z -S3 S4C3 C4C3 Z 

( 12 ) 
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Consider a unit vector along the X body axis. This will yield the following 

x r Y r Z I 


IC 


IC 


IC 


components. 


= C3CL 

(13) 

= C3SL 

(14) 

= -S3 

(15) 


Define the pitch angle of the vehicle with respect to the Xj, Y^, frame to 

be a . Define the bearing angle of the vehicle with respect to the reference 
bearing line to be /3. It follows that 

" ZlC 

tan a = — 

\[~X? + Y 2 
'IC IC 

where X , Y and Z are given in equations (13), (14), and (15), 

1L 1L IL 

respectively. Therefore, 

S3 

tan a = — ■-■■■ — ■- 

n/ C3 2 CL 2 + C3 2 SL 2 

or 

tan a = tan 3 


and 


angle a = angle 3 


(16) 
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Although this result is no great surprise, it is significant that the pitch angle 
can be determined by direct reading of gimbal angle A3. 

The bearing angle is given by 


tan (3 



tan L 


or 


angle j3 = angle L 


From equations (8) and (9) one finds that 


angle L = tan 


-i 


C1S2C4 + S1S4 


C1C2C3 - C1S2S4S3 + S1C4S3 


(17) 


so that the bearing angle is defined in terms of the gimbal angle readings 
Al, A2, A3, and A4. 

Another angle of interest is the angle from the vehicle heading reference 
to the point location of the LRV projected onto the X^., plane. This angle 

is illustrated in Figure 4. It is desired to know the angle that the vehicle must 
head to return to the LEM location. We define this angle in Figure 5. The 
angle $ , defined as the counterclockwise rotation from the X , projected 

heading of the vehicle, gives the desired angle. For the general case. 


j3 + <p - e = 180 degrees 


or 


<f)= 180 + 0-/3 


(18) 
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where 


e 


tan 1 



( 19 ) 


With the information given in equations (16), (17), and (18), navigation 
requirements can be satisfied. A block diagram of the navigation system is 
shown in Figure 6. For convenience the displayed quantities are defined. 


S 

s 

R 

P 

<P 


- Total distance traveled by the LRV 

- Velocity of LRV 

- Radial distance traveled by the vehicle 

- Bearing angle of the vehicle with respect to the reference 
bearing 

- Heading angle of the vehicle back to LEM 


Xj, Yj, Zj - Distance traveled by the vehicle with respect to the moon- 
fixed inertial system whose origin rotates with the moon. 


SIMULATION CONSIDERATIONS 


To determine the range of required quality of the gyros and acceptable 
system misalignments, error analysis for the navigation system must be made. 
The details of how a digital simulation of the navigation system can be used to 
perform this error analysis are as follows. 


Outline of Simulation Study 

We first consider the effects of gyro drift errors. These errors, when 
multiplied by time, are reflected in the gimbal angle readings of Al, A2, A3, 
and A4. The following method is used to assess the effects of gyro drift rates. 
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FROM GYROS 



Figure 6. Block diagram of proposed navigation system 













Step 1 . > Assume a velocity profile for the vehicle. This velocity 
is denoted by S. 

Step 2. Assume a set of sequential angles, in particular the angle 
sequence given in Figure 3, corresponding to the velocity profile in Step 1. 

Step 3. From the rotation sequence in Step 2, develop the expressions 
for determining Al, A2, A3, and A4 as a function of the angles in Step 2. 

Step 4. Write a digital simulation program for. the navigation system 
given in Figure 6. The inputs to the program will be S and Al, A2, A3, 
and A4. 


Step 5. Run the simulation for a specified time, T The following 

information will be printed out every AT seconds: time, X^, Y^, Z^, /3, 4 > > 

and R. These printouts will establish a reference set of data to be compared 
with future runs which include drift rate errors. 


Step 6. Introduce a set of. known drift rates into Al, A2, A3, and A4 
as given in Step 3. 

Step 7. Run the simulation for the same velocity profile as given in 
Step 1. Print out the same data at the same AT as produced in Step 5. 

Step 8. Compare the results of Step 7 with Step 5 to determine the effect 
of errors introduced by the drift rates. These errors should be compared with 
the required system accuracy specifications as previously given. 

Step 9. Repeat Steps 6, 7, and 8 for as many combinations of drift 
rates as desired. 

To assess the effects of misalignments, a similar approach to that 
outlined in Steps 1 through 9 can be employed. There are two major sources 
of misalignment errors. First, the two gyros, may be physically misaligned 
with respect to the vehicle frame. The misalignment of the upper gyro is 
independent of the lower gyro. Second, on any given sortie, the vehicle may 
not be aligned precisely as it was on previous sorties. This means that an 
error will exist with respect to locating previously established sites. 
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Hardware misalignments can be held to within 30 arc seconds by 
careful alignment procedures. Since these errors are small compared to 
vehicle misalignments, they will be neglected. 

Vehicle misalignment errors can easily be evaluated by specifying the 
sequential angles L, A3, and A4 corresponding to the desired misalignment 
to be evaluated. These misalignment angles are introduced in Step 6 of the 
error evaluation procedure. 


Background Analysis for Simulation Inputs 

This section concerns the problems of selecting a velocity profile S 
for Step i, selecting a body rotation sequence (as a function of time) for Step 2, 
and developing the expressions for determining the angles Ai, A2, A3, and 
A4 as a function of the sequential angles. First, the relationships between the 
angles will be developed. 

The rotation sequence selected is given in equation (i) and shown in 

Figure 3. The first rotation is around the Z axis by an angle L, the 

13 

second rotation is around the Y axis by an angle A3, and the third rotation 

13 

is around the X axis by an angle A4. Note that the sequential angles A3 and 

A4 are the same angles A3 and A4 as read by the lower gyro gimbals in 
Figure 2. It has been shown that if angles Ai, A2, A3, and A4 are given, 
then the sequential angles L, A3, and A4 can be determined. Actually the 
only angle to be determined is angle L, and it may be determined from 
equations (8) and (9). Now to make the relationship between the sequential 
angles and the gimbal angles bilateral, expressions for AI and A2 must be 
determined as a function of angles L, A3, and A4. 

The transformation matrix between the inertial axis and the body axes 
can be written as 


T 


C2C1 C2S1SU + S2CU 
-S2C1 -S2S1SU + C2CU 
Si -CiSU 


-C2S1CU + S2SU 
S2S1CU+ C2SU 
CiCU 


( 20 ) 
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or 



C3CL 

C3SL 

-S3 

T = 

-C4SL+ S4S3CL 

C4CL + S4S3SL 

S4C3 


S4SL + C4S3CL 

-S4CL+ C4S3SL 

C4C3 


Therefore, equating equation (20) to (21) and solving for Ai and A2 yields 

Ai = sin -1 [S4SL + C4S3CL] (22) 


and 


A2 


tan -1 


[C4SL - S4S3CL] 
C3CL 


(23) 


Thus, with knowledge of the sequential angles L, A3, and A4, the gimbal 
angles AI, A2, A3, and A4 are known and this completes the bilateral 
relationship. 


Vehicle Velocity and Position Profile for Error Analysis 

The vehicle velocity and position profile can be arbitrarily selected. 
However, this method may lead to unrealistic vehicle position. For example, 
the vehicle may be traveling upside down. To avoid this, a profile will be 
selected so that the position and velocity of the vehicle are approximately 
known as a function of time. The profile selected is shown in Figure 7 and 
is divided into five segments as indicated. The following information describes 
each segment. 

Segment i. Over segment 1 the vehicle has a linear velocity given by 


S = 30t km/hr (0 t 0. 3 hr) 
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Figure 7. Selected vehicle profile for error analysis. 



The heading of the vehicle is at a constant 30 degrees so that 

angle L = +30 degrees 
angle A3 = 0 degree 
angle A4 = 0 degree 

The total distance traveled is 

S = i. 35 km 


At t = 0. 3 hr, the velocity of the vehicle will be 9 km/hr and the coordinates 
of point A will be 


Xj = 1.17 km, Yj = 0,675 km, = 0 km 

The gimbal angles A1 and A2 are determined from equations (22) and (23) , 
respectively. 


Ai = 0 degree, A2 = 30 degrees 


Segment 2 . Over segment 2 the vehicle has a constant velocity of 


S = 9 km/hr (0. 3 < t < 0. 796 hr) 


The vehicle has a constant 0-degree heading, but a sinusoidal pitch profile 
as shown in Figure 8. The distance traveled along the sinusoidal path is 
given by 

12 o' 

S = — / v 1 + K| M 2 cos 2 MX d(MX) (22) 

o 
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Figure 8. Pitch profile for segment 2 of vehicle path. 

where K* = 0. 2 km; M = 5.23. Equation (22) can be put into the elliptic 
integral form 

o’ i 

S = A J n/ i - K 2 sin 2 !/) d ip 

o 

where 


K 




1 + kf 




9 


A 


12 ^ 1 + Ki M 2 
M 


The approximate evaluation of this elliptic integral of the second kind gives 
S = 4. 464 km 


The sequential angles are given by 


angle L 
angle A4 
angle A3 


0 degree 
0 degree 

46. 3 degrees cos 


0.744 


t» 


) 
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The gimbal angles A1 and A2 are 
Ai = 46.3 degrees cos 
A2 = 0 

where 

(t* = t - 0. 3) 

The Xj, Yj, Zj coordinates at point B are approximately 

Xj = 4.77 km, Y J = 0.675 km, Zj = 0 km . 

Segment 3 . The velocity of the vehicle is constant over this segment 
and given by 

S = 9 km/hr (0. 796 < t < 0. 996) 

The heading of the vehicle is a constant -90 degrees. The distance traveled 
is 

S = 1. 8 km 

The vehicle has a sinusoidal roll motion of 
30 degrees sin (1007rt") 

t" = t- 0.796 (0.796 ^ t< 0.996) 
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The sequential angles are given by 


angle L = -90 degrees 

angle A3 = 0 degree 

angle A4 = 30 degrees sin (lOOirt") 
The gimbal angles A1 and A2 are 


Ai = -30 degrees sin (100?rt n ) 

A2 = -90 degrees 

The Xj, Yj, Zj coordinates at point C are 

Xj = 4.77 km, = -1.125 km, Z^ = 0 km 

Segment 4 . The vehicle velocity is constant over this segment and 
given by 


S = 5.01 km/hr (0. 996 t < 2. 39) 

The pitch and roll angles of the vehicle are zero and the yaw profile is shown 
in Figure 9. The length of the path can be evaluated as an elliptic integral of 
the second kind and yields 


S = 7. 0 km 

The sequential angles are given by 

angle L = 180 degrees - 57. 5 degrees cos (31. 506t m ) (t m = t - 0.996). 
angle A3 = 0 degree 
angle A4 = 0 degree 
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Figure 9. Yaw profile for segment 4 of vehicle path. 

The gimbal angles Ai and A2 are 

A1 = 0 degree 

A2 = 180 degrees - 57. 5 degrees cos (3i.506t ,M ) 

The Xj, Yj, Zj coordinates at point D are approximately 

Xj = 0 km, Yj = -1. 125 km, Zj = 0 km 
* 

Segment 5 . The vehicle velocity is constant over this segment and 
given by 

S = 1.85 km/hr (2.39 <t< 3.0) 


23 




The sequential angles are 


angle L = 90 degrees 
angle A3 = 0 degree 
angle A4 = 0 degree 

The gimbal angles Ai and A2 are 

Ai = 0 degree 
A2 = 90 degrees 

The Xj, Yj, Zj coordinates at point E are approximately all zero. 

The answers for all segments above are only approximate. The actual 
values with zero system errors are given by the computer simulation results. 


A DIGITAL SIMULATION OF THE NAVIGATION SYSTEM 


A block diagram of the navigation system is shown in Figure 6. This 
block diagram is used as the starting point for writing the simulation program. 
It should be pointed out that Figure 6 shows the complete set of calculations 
required for this navigation system. On the other hand, the digital simulation 
will require additional calculations. In particular the angles AI, A2, A3, and 
A4 are direct inputs to the actual system, but these same quantities must be 
calculated as a function of time in the digital simulation. Also, a capability 
must be included in the digital program for simulating system errors. 

A flow chart of the simulation program written for this navigation 
system is shown in the Appendix. The actual program was written in Fortran 
and the quantities given in the program listing (Appendix) were modified to 
be compatible with Fortran requirements. 
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SIMULATION RESULTS 


Using the surface profile shown in Figure 7, a simulation run was first 
made without system errors to establish a set of reference data. The results 
are in good agreement with the analytical approximations given previously 
and therefore validate the computer simulation. 

Graphs showing the reference data for X^, Y^, Z^, fi, <fi, and R 

versus time are given in Figures 10 through 15, respectively. For convenience, 
the definitions of these quantities are given below. 


Xj, Yj, Zj - Distances traveled by the vehicle with respect to the 
moon-fixed inertial coordinate system. The origin of 
the moon-fixed coordinate system is located near the 
LEM and moves with the lunar surface. 


P 

<P 

R 


The angle between the reference direction (X^ axis) and 
the direction of travel of the vehicle. 

The angle between the direction of travel of the vehicle 
and the shelter site. 

The straight line distance between the shelter site and 
the LRV. 


All subsequent computer simulation runs with gyro drift rates and vehicle 
misalignments are compared to the reference data. 

Table 1 gives a summary of the simulated system errors. Various 
cases for gyro drift rates and vehicle misalignments were simulated. These 
cases are representative of the expected range of errors but do not include 
all possible combinations in this range. To cover the entire spectrum would 
require an extensive number of computer runs. However, from the cases 
studied one can establish trends on the accuracy of the navigation system as 
a function of system errors. 
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Beta versus time for reference trajectory. 
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Figure 15. Range versus time for reference trajectory. 
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TABLE 1. SYSTEM ERRORS FOR SIMULATION STUDY 
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Figure 16 shows a plot of the error in returning to the LEM site as 
a function of gyro drift rates. Drift rates in gimbal angles Al, A2, A3, and 
A4 are of equal magnitude and sign. The end point error is practically the 
same for either positive or negative drift rate polarity. The error in returning 
to the shelter site varies linearly with the gyro drift rates. The diagram shows 
that the vehicle will return to within a 500-meter spherical range (the maximum 
allowable error) for drift rates up to 1. 5 degrees/hour. However, if only 
in-plane errors are considered, the vehicle will return within a 500-meter 
circle for drift rates up to 4. 8 degrees/hour. For all cases, the major 
contributing factor to the end point error is the Z^ component. In making 

any closed trajectory, the error accumulation in the axis is monotonic. 

On the other hand, the accumulative errors in the and axes have a 

tendency to balance out because the trajectory returns to the starting point, 
forming a closed path. 

From a physical standpoint, the vehicle is restricted to the surface of 
the moon. If the endpoint errors in and are small, the vehicle will 

be physically close to the starting point regardless of Z^ errors. However, if 

one is attempting to determine the altitude of the vehicle, errors in Z^ become 
serious. 

In Table i, cases i through 12 are for gyro drift rates only. Vehicle 
misalignments are zero for these cases. Figures 17 through 28 show the 
error in the angle back to the LEM site for the respective cases. The nature 
of this error angle should be explained. 

For the first simulation run, the vehicle follows the reference 
trajectory with a given velocity profile without any system errors. The angle 
<p from the heading of the vehicle to the LEM site is calculated as a function 
of time. A second simulation run is made with gyro drift rates (multiplied 
by time) added to the gimbal angle readings. In this case the vehicle physically 
follows the same reference trajectory of the first simulation run. In other 
words, the vehicle pilot can visually navigate, without the aid of the navigation 
system, along the tracks made by the reference trajectory. Consider that 
the vehicle is within visual range of the LEM and located at point A on the 
reference trajectory shown in Figure 29. The pilot can observe that he must 
turn the vehicle through an angle <p ^ to head directly to the LEM. However, 

the information from the navigation system may locate the vehicle at point B 
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Figure 17. Bearing angle error for drift rates of A1 = +1, A2 = +1, 
A3 = +1, A4 = +1 deg/hour. 



Figure 18. Bearing angle error for drift rates of A1 = +2, A2 = +2, 

A3 = + 2 , A4 = +2 deg/hour. 
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<r drift rates of Ai = +3, A2 = +3, A3 = +3, A4 = +3 deg/hour. 



END OF SEGMENT 
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Figure 21. Bearing angle error for drift rates of Ai = -2, A2 = -2, A3 = -2, A4 = -2 deg/hour. 


END OF SEGMENT 



Figure 22. Bearing angle error for drift rates of A1 = -3, A2 = -3, A3 = -3, A4 = -3 deg/hour. 



Figure 23. Bearing angle error for drift rates of A1 = +1, A2 = -1, 
A3 = +1, A4 = -1 deg/hour. 



Figure 24. Bearing angle error for drift rates of Ai = -i, A2 = +1, 

A3 = -1, A4 = +1 deg/hour. 






Figure 27. Bearing angle error for drift rates of Ai = -2, A2 = + 2 , 
A3 = -2, A4 = +2 deg/hour. 



Figure 28. Bearing angle error for drift rates of Al = -2, A2 = - 2 , 

A3 = + 2 , A4 = +2 deg/hour. 
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and give a bearing angle to the LEM of <f> . The error in the bearing angle 

is defined as 0 E ~ if the driver turns the vehicle through an angle 

4> with the vehicle located at point A, the error in the bearing angle is <p . 

J3 E 

The angle <j> versus time is plotted in Figures 17 through 28. As the 

vehicle approaches the LEM site, the bearing angle error always increases 
in magnitude. The reason for this increase is easily understood by referring 
to Figure 29. When the range distance becomes small (less than 300 meters) , 

4> can be considerably greater in magnitude than <fi . . Although the resulting 

error, is large, it is not significant because the vehicle driver can 

easily observe the location of the LEM within the 300-meter range. 

The second area of interest in the simulation study is the effects of 
vehicle misalignments. The ability to revisit an established site depends upon 
how well the vehicle can be realigned to the original reference direction. The 
effects of misalignments can easily be evaluated by intentionally offsetting the 
sequential angles A3, A4, and L. Cases 13 through 23 in Table 1 give the 
misalignments that were simulated. The misalignment errors were evaluated 
by the instantaneous radial error expressed by 

* E T (t > - 4[ X I (t . ) - X IE ( V] 2+ [ Y I ( V - Y IE<V ] 2+ [ Z I (t a> - Z IE (t a>] ' 

t=t 

a 
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where X T (t ) , Y (t ) , Z (t ) are the coordinates along the reference path at 

i 3. 1 cl 1 cl 

t = t , and X TT , (t ) , Y t _, (t ) , Z T _, (t ) are the coordinates at t = t as given 
a IE a IE a IE a a 

by the navigation system when misalignment errors are added. Figures 30 

through 40 show the variations in AR^. as a function of time. AZ is the most 

significant component of AR^ and represents the vehicle altitude error. If 

altitude errors are neglected (this is reasonable since the vehicle must remain 
on the lunar surface) , the important error is the in-plane error given by 


AR Tp (t) 


t=t 


^[w- x *“•>] !+ [ 


w - w 


a 


This error, as a function of time, is illustrated in Figures 30 through 40. 



Figure 30. Radial path error versus time for vehicle misalignments 
of A3 = + 1 deg, A4 = +1 deg, L = +1 deg. 

Nominal in-plane errors are noted from the diagrams to be generally 
less than 200 meters. The system requirement is that the vehicle must have 
navigation capability to revisit an established site. The ability to revisit a 
site depends upon land identification marks in the site area. Site identification 
can be established by flag markers or low level signal radio transmitters. In 
any case the navigation system must be capable of returning the vehicle to 
within identification range. A navigation system with a site return accuracy 
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Figure 3i. Radial path error versus time for vehicle misalignments 
of A3 = -1 deg, A4 = -1 deg, L = -1 deg. 

of 200 to 300 meters should be adequate for relocating any given point within 
the maximum 5-km radial distance from the LEM. The preliminary simulation 
results indicate that misalignments should be less than 5 degrees to return to 
within 300 meters of any established site on the path profile. 


CONCLUSIONS 


The basic hardware required to implement the navigation scheme for 
lunar exploration is two free gyros mounted to the frame of the vehicle. A 
particular lunar trajectory was selected to evaluate the accuracy of the system. 
The time varying sequential angles were used as inputs for a digital simulation 
to study gyro drift rates and vehicle misalignments. The results of the simu- 
lation study indicate that system requirements can be satisfied with this navi- 
gation scheme . To meet the requirement of ±3 degrees pointing accuracy to 
the LEM site, the gyro drift rates should be in the range of 1 degree/hour. 

The capability to revisit any point (within 300 meters of accuracy) requires 
vehicle realignment accuracy of 5 degrees or less. 

This navigation scheme represents only one of many possible candidates. 
Particular hardware configurations that should be analyzed include a directional 
and vertical gyro system, a sun aspect sensor and vertical gyro system, and a 
directional gyro and inclinometer system. The navigation equations for each 
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Figure 32. Radial path error versus time for vehicle misalignments 
of A3 = +3 deg, A4 = +3 deg, L = +3 deg. 




END OF SEGMENT 



Figure 33. Radial path error versus time for vehicle misalignments 
of A3 = -3 deg, A4 = -3 deg, L = -3 deg. 
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Figure 34. Radial path error versus time for vehicle misalignments 
of A3 = 1 deg, A4 = 1 deg, L = 1 deg and drift rates of i deg/hour. 
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Figure 35. Radial path error versus time for vehicle misalignments 
of A3 = +3 deg, A4 = +3 deg, L = 3 deg and drift rates of +1 deg/hour. 





Figure 36. Radial path error versus time for vehicle misalignments 
of A3 = i deg, A4 = i deg, L = 1 deg and drift rates of -1 deg/hour. 
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Figure 37. Radial path error versus time for vehicle misalignments 
of A3 = +3 deg, A4 = +3 deg, L = +3 deg and drift rates of -1 deg/hour. 
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Figure 38. Radial path error versus time for vehicle misalignments 
of A3 = +1 deg, A4 = -2 deg, L = 3 deg and drift rates of +i deg/hour. 
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Figure 39. Radial path error versus time for vehicle misalignments 
of A3 = +5 deg, A4 = +5 deg, L = +1 deg. 





Figure 40. Radial path error versus time for vehicle misalignments 
of A3 = + 4 deg, A4 = +10 deg, L = 1 deg. 





of these systems should be developed and related directly to measurable angles 
from the hardware. Simulation studies should be made to determine the system 
accuracy. Finally, a system selection should be made which represents the 
best compromise between cost, weight, reliability, minimum component develop- 
ment, and navigation accuracy. 
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APPENDIX 


Figure A-l gives the flow chart of the navigation system simulation. 
The program listing is given in Figure A-2. 



Figure A-l. Flow chart of navigation system simulation. 
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C LUNAR RBVER VERS I BN TWO 

DIMENSI8N A(4)/Al(4)/A2(4) 

DATA T/O*/ 

Data xi/q*/ 
data yi/o#/ 

Data zi/o«/ 

Data DELTA/ #00390625/ 

Data s/o*/ 

Data r/o#/ 

Data vi#v2/o##o#/ 

Data fee/o./ 

C 6D8NETER 

C BEARING and inclination 

call headng 

WRITE. (5/1) T/S/V2/ZI/XI/YI/R/FEE/B 

call velqc<T/vi) 

2 call velqcc (t+delta)/V2) 
call gimbaut/A) 
call errqr (t>a) 
call bearng(A/B) 

DSs( (V l+V2)/2. )*DELTA 
S*S+DS 

fc coordinate displacements 

ZI*ZI+DS*S1N(A<3) ) 

HD«DS*CSS(A(3) ) 

XI s Y I +MD*CBS ( B ) 

Vl«Yi+HD*SlN(B) 
c Range 

R*sqrt < X I *X i +Y I * Y I +Z I *Z I ) 

C DIRECTION T8 hbme base 

Call angle(yi/XI/theta) 

FEE«3*1A1592653589793+THETA'B 

t«t+delta 

C 8UTPUT TB GUAGES 

WRITE (5/1) T/S/V2/ZI/XI/YI/R/FEE/B 

IF(T*GT»3*10) CALL EXIT 

Vl*V2 

DS 6 K-l/4/1 
6 A1(K)»A2(K) 

GB TB 2 

1 FBRMAT(9(F12#S/2X) ) 

END 


Figure A-2. Simulation program listing. 
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SUBROUTINE VELOC ( T> V ) 

DIMENSION TR(5) 

Oat a TR ( l ) t TR < 2 ) * TR (3 ) , Tr ( 4 ) , TR ( 5 ) / * 3/ * 792* #994* 2 *39, 3 . 0 / 
If d # le*tr( 1 ) ) Ge to 1 

IF ( T »LE»TR (2) ) G9 T9 2 
IF ( T»LE *TR ( 3 ) ) GB T9 3 
IF (T#LE*TR(4) ) GB T9 4 
IF ( T #LE*TR ( 5 ) ) GB T9 5 
V« 0* 

Return 

1 V*30*#T 

Return 

2 V®9 *6 
Return 

3 V*9»0 
Return 

4 V?4# 60983 

Return 

5 V* l » 8851 

Return 

end 


SUBROUTINE ERR0R(T*A) 

dimension a(4) 

Return _ 
end 


SUBR8UTINE GlMBAL(T*A) 

dimension aui 

CALL DR I VER ( T * A * AL ) _ 

Sai«(SIN(A( 4) )#SlN(AU+CeS(A(4) ) *SI N( A ( 3 ) ) *C0S ( AL ) ) 

CA 1 *SQR.Ti I • *SA1 *SA1 L 

Call anglE(Sai*cai*a(i ) ) 

SA2*CBS(A(4) ) *S I N ( AU ) "S I N ( A ( 4 ) ) *SIN < A ( 3 ) )*C0S(Ak> 

CA2®C0S(A(3M*C0S(AU 

CALL ANGLE ( SA2* CA2* A ( 2 ) ) 

Return 

end 


Figure A-2. (Continued) 
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SUBR9UTJNE DRIVER(T/A,AL) 

D 1 MENS I ©N TR ( 5 ) . A ( 4 ) 

Data TR(1)/TR(2)/TR(3)/TM4)>TR<5)/.3/ .792/.994*2.39/3.0/ 
DATA FAC/O. 017453292519943/ 

Data PI/3.141592653589793/ 

Data W2/38. 2085592993/ 

Data W3/314.1592653589793/ 

IF(T.UE.TRU) ) Ge T0 1 
IF ( T*l_£*TR ( 2 ) ) G9 T9 2 
IF(T.UE*TR(3) ) 39 T9 3 
IF ( T tLE.TR ( 4 ) ) G8 T9 4 
IF ( T *LE»TR (5 ) ) G8 T9 5 

return 

1 AL*. 5235988 
A ( 3 ) *0 * 

A ( 4 ) *0* 

Return 

2 AL*Qt 

A ( 4 ) »0 • 

A(3)sFAC#46*29*CeS(W?*(T-TR(l) ) ) 

Return 

3 Al = n. 5707963 
A ( 3 ) * 0 # 

A(4)».5235988*SIN(W3#(T-TR(2) ) ) 

Return 

4 Au s Pl^57.4*FAC*CeS(31.5Q59435164* (T*TR(4)>> 

A ( 3 ) *0 * 

A ( 4 ) *0* 

Return 

5 AU-1, 5707963 
A ( 3 ) «0« 

A ( 4 ) «0* 

RETURN 

end 


Figure A-2. (Continued) 
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Figure A-2. (Concluded) 
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